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We investigate the decoherence process for a quantum register composed of N qubits coupled to 
an environment. We consider an environment composed of one common phonon bath and several 
electronic baths. This environment is relevant to the implementation of a charge based solid-state 
quantum computer. We explicitly compute the time evolution of all off-diagonal terms of the 
register's reduced density matrix. We find that in realistic configurations, "superdecoherence" and 
"decoherence free subspaces" do not exist for an iV-qubit system. This means that all off-diagonal 
terms decay, but not faster than e~ q ^ N , where q(t) is of the same order as the decay function of a 
single qubit. 
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I. INTRODUCTION 

A typical quantum computer would consist of a large 
number (iV) of two-level quantum systems, coined qubits, 
where the level splitting of each qubit and the interac- 
tion between pairs of qubits is adjustable. Quantum op- 
erations are then obtained by varying these parameters 
along a scheme defined by a quantum algorithm. The 
physical system composed by the N two level systems is 
our quantum register. In the ideal case, when the quan- 
tum register is isolated, the time evolution of an arbitrary 
initial state ipi of the register is unitary. Such an ideal 
quantum computer could be used to solve some problems 
more efficiently than classical computers. An important 
example is Shor's quantum algorithm to factor an integer 
with n digits in a time growing polynomially with n in- 
stead of exponentially when using a classical computer. 1 

However, any realistic quantum computer is coupled 
in some way to an external environment, which leads 
to decoherence and dissipation. The quantum register 
becomes entangled with the environment, so that its ef- 
fective evolution is not unitary anymore. There can also 
be energy transfers between the register and the environ- 
mental bath, which lead to dissipation and decoherence. 
However, dissipation is not a requirement for decoherence 
to occur. 2 

The quantum decoherence process is elegantly ex- 
pressed in the framework of the reduced density matrix 
of the quantum register. When no coupling to the en- 
vironment is present, the reduced density matrix simply 
follows a Heisenberg-type evolution. As soon as the cou- 
pling to the environment is introduced, the off-diagonal 
terms of the reduced density matrix of the register decay 
with respect to time. This is often referred to as phase 
damping. In the simplest case of a single two level system 



connected to an environment, the off-diagonal elements 
of the reduced density matrix decay exponentially in time 
as ~ e _Q, w, where t is the time and the function q(t) de- 
pends on the strength of the coupling to the environment. 
In the context of quantum information processing, such 
a decoherence event can be expressed as a quantum er- 
ror. Following a pioneering work by Shor 3 , it was shown 
that it is possible to encode each qubit using a minimum 
of five qubits (the five qubit code) in conjunction with 
quantum error correction algorithm in order to "repair" 
a faulty qubit. 4 This would enable an accurate quantum 
computation as long as the error rate is small. The draw- 
back of all quantum error correction schemes is that at 
least five times as many qubits are necessary for the same 
operation than in the ideal case. 

In the case where there are N two level systems the 
situation is potentially much worst since the decoherence 
of the register cannot be simply expressed as a super- 
position of single qubit decoherence. Indeed, Palma et 
al. argued that the decay of the most off-diagonal ele- 
ments goes like ~ q-iw-N (superdecoherence), when all 
qubits are imbedded in a single bath. 5 Such a dependence 
would jeopardize the use of quantum error correction al- 
gorithms as soon as q(t)N approaches 1, since the error 
rate would simply become too large. A potential res- 
cue to this problem was proposed with the existence of 
decoherence free subspaces. 6 

In this article we investigate the decoherence of N 
qubits in the context of a realistic collection of solid state 
two level systems (our qubits) imbedded in a semicon- 
ducting environment. More specifically, we consider the 
case where the two-level system is a single electron in 
a double quantum dot patterned in a two dimensional 
electron gas confined in a GaAs/AlGaAs heterostructure. 
These qubits are coupled to a common phonon bath and 
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to additional independent electronic baths, representing 
the metallic leads which allow the control and opera- 
tion of the qubits. Experimentally, a decoherence time 
of about 1 ns was recently measured in a single double 
quantum dot. 7 For this system, we find that the deco- 
herence, i.e., the decay of the most off-diagonal elements 
of the JV-qubit reduced density matrix, goes like e~ q ^ N , 
a much slower decoherence rate that previously thought, 
but that decoherence free subspaces do not exist. 

To obtain these results, we consider a model based on 
a scaled version of the experimental double quantum dot 
as seen in Fig. 1. Tanamoto showed that such a system 
can perform all the operations necessary for a quantum 
computer. 8 Many other groups have also used similar 
coupled quantum dots geometries as model system for 
a qubit. 9-18 While we use this particular system for our 
model, our results are in fact more general and remain 
qualitatively similar for different physical realizations. 
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FIG. 1: Schematic representation of a solid state quantum 
computer coupled to a common phonon bath and several in- 
dependent electronic baths. 

The rest of this article is organized as follows. In 
Section 11 we introduce the Hamiltonian, similar to a 
spin-boson model, and describe the ./V-qubit register, the 
environment and the coupling between the register and 
the environment. In Section III we give an exact formal 
expression for the reduced density matrix of the regis- 
ter. The time evolution of the reduced density matrix is 
expressed in terms of the influence functional, which is 
computed in Section IV. In these sections the expressions 
are fairly general and do not depend on the exact model 
considered. In order to gain insight into a physical sys- 
tem we consider the system shown in Fig. 1 in the next 
sections. Sections V and VI are devoted to the special 
case of acoustic phonons coupled to charge-qubits. The 
coupling to electronic baths is analyzed in section VII. 
In Section VIII we evaluate the decoherence rate for the 
coupling to a single phonon bath, where we show that 
the decoherence function scales as e~~ q ^ N when increas- 
ing the number, N of qubits. We give physical estimates 



for piezo and deformation phonons in Section IX and 
analyze our results in the dynamical case in Section X, 
where we introduce quantum operations on the register 
and evaluate the decoherence process. Finally, a short 
summary and conclusions are given in the last section 
XL 



II. THE HAMILTONIAN 

We use a Hamiltonian for our model qubit, which de- 
scribes the tunneling of a single electron tunnelling be- 
tween two adjacent quantum dots. The electronic state of 
the dots can be controlled by adjusting the gate voltage, 
which allows either localization of the electron or reso- 
nant tunnelling between the dots. The complete physical 
localization of the electron in a given dot is denoted by 
the vector (1,0) T while localization in the other dot cor- 
responds to (0,1) T . At low enough temperatures, only 
the combinations of these two states need to be taken into 
account. Each qubit is described by the Hilbert space C 2 
and the canonical basis is denoted by | + 1) = (1, 0) T and 
| — 1) = (0, 1) T . The single qubit Hamiltonian reads (we 
write the Schrodinger equation as iip — Hip, hence the 
units of H are s _1 ) 



H = -A t cr x 



£tcr z 



(1) 



where a x and a z denote the Pauli's matrices, u z \ ± 1) = 
±| ± 1) and <r x \ ± 1) = | =F 1) . Typically, the tunnelling 
matrix element is A ~ 1 — 100GHz and the bias e is 
adjusted with the gate voltage. Both quantities need to 
be dynamically controlled for the operation of the qubit 
in a quantum computer. 8 

For the ./V-qubit register, we write |^) (where I S 
{-1,+1} W ) for the vector of {C 2 )® N defined as \l) := 
\h) (8 • • • <8 \In)- The Hamiltonian of the register is de- 
composed as H r (t) = £(t) + A(t), where £(£) denotes 
the diagonal part of H r (t) with respect to the basis \l), 
that is £(i)|Z) = e(t,l)\l). 

For the total system, i.e., the register composed of 
N qubits plus environment, we consider the following 
Hamiltonian 



H{t) =H r {t) + H e + H r 



H' 



(2) 



<M>L + d>lb 



k°k- 



where (f>k are operators acting on the register's Hilbert 
space. We assume that they are diagonal with respect to 
the basis \l). The field operators 6k and are bosonic 
operators, that is 



bir> , b 



(3) 



and all other commutators are zero. The Hamiltonian 
given in Eq. (2) is very general and is frequently used 
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to study open quantum systems. It leads directly to the 
spin-boson model, which describes many types of envi- 
ronments with extended degrees of freedom (phonons, 
electrons, magnons...) and can be obtained from micro- 
scopic models. 5 ' 19,20 It is however inappropriate for lo- 
calized environments, such as a bath of nuclear spins. 21 
We decompose H(t) as H(t) = H (t) + A(t) with 



H {t) = H{t) + H e + H re . 



(4) 



Note that (l\H (t)\m) = 8 lm (e(t,m) + H e + H re (m)) 
where 



k 

0k (m) = (m|0k|m) . 



(5) 



As a consequence, we find that the evolution operator 
associated to Ho(s) is given by 

(l\U {t,0)\m) = 5 lme -itic{s,™)d8 e -H(H°+H™ {m )) (6) 



III. EXACT FORMAL EXPRESSION FOR THE 
REDUCED DENSITY MATRIX 

We now give an exact formal expression for the reduced 
density matrix of the register by expanding the evolution 
operator of the total system with respect to A(i), the off- 
diagonal terms of H r (t). This corresponds to the method 
of Leggett et al studied at length for the spin-boson model 
in Ref. 19 . 

The density matrix of the total system is given by 
p(t) = U(t,0)p U(t,0)i where po is the density ma- 
trix at time t = and U(t, 0) is the evolution opera- 
tor associated to the total Hamiltonian H(s). We de- 
fine the interaction picture with respect to Ho(s) as 
p{t) = Uo(t,0y p(t)Uo(t,0). Hence, the Heisenberg equa- 
tion reads ip(t) = L(t)p(t), where L(t) is the Liou- 
ville operator Lit) A = [A(t),A]. Therefore we have 
pit) = [/ (t,0)(T{exp[-i f*Lis)ds}}p Q )Uoit,Oy. 

The density matrix of the register is defined by trac- 
ing out all the environment degrees of freedom p r it) = 
Tr e [p(t)]. So, we need to compute the trace over the de- 
grees of freedom of the environment of terms of the form 



(l\Uoit, 0)A(^_ j ) • • • A(t{)p A(ti) • • • A(tj)tf (t, 0)», 

(7) 



4> 



with t[ < t 2 ■ ■ ■ < t p _ j and t{ < t l 2 ■ ■ ■ < t\. This term 
reduces to a sum over the all pair of maps {C} , C}) defined 



on [0, t] with values in { — 1, +1}^, constant by step, and 
with (^it~) = I and CH* _ ) = m - More precisely, C} is 
constant on each interval [tj, where < r < p — j, 
with 4 = and t\_ j+1 = t. Moreover, C T (*J>-j) = 1 0n 
the other hand, is constant on each interval [4)4+1 1 
where < s < j, with 4 = and tjj +1 = t. Moreover, 
C^tj) = m. We write C r T for C T (4), and Cj for CH4)- 
With these notations, the term above becomes 



E II ri(cJ|A(4)lc T -i)(d-i|A(4)ici) 

(fi,£T) r=l s=l 

x(Cl_ j \Uoit,tl_ j )\Cl_ j ) ■ ■ ■ (CoVo(4,0)|Co T ) 

x(Co T ipoiCo i )<Co i i^o(o,4)iCo i }---<cii^(ti^)icj) . 

(8) 

We now assume the register and the environment to be 
decoupled at time t = 0, hence p — Pq <8> p%- Moreover 
the environment is assumed to be initially at thermal 
equilibrium which means that p% = e-WZ' 1 , where 
= h/K B T, hence 0T ~ 7.64 • 10~ 12 sK. 

Let 7 be the path in the complex plane defined as 
j(s) = s for s e [0,t], 7(s) = 2t - s for s &]t,2t], and 
7(s) = -i(s - 2t) for s e]2i,2i + /?]. Define £(s) as 
C(«) = C T (7(*)) for s G [0,t], and C(«) = C'M*)) for 
s 2i], and £(s) = for s e]2t, 2t + /?]. For A > 0, we 
define 

C/a( S ) = t| C x P -i^ S 7 ( S ')(ff e + Ai/ re (C( S ')))^' ]• 

(9) 

Moreover, we define the influence functional as 



Z x [() = —Tr e [U x i2t + p)}. 



(10) 



Finally, if p and j are integers with < j < p, we 
denote by O w - the set of all pairs (£7i, SI2) of subsets 
of {1, •••,£>} such that f^i = {pi,---,Pj} and = 
{qi, ■ ■ ■ , q P -j} with pi < ■ ■ ■ < Pj, and qi < ■ ■ ■ < q p -j, 
and such that Oi U Cl 2 = {1, • • • ,p} (hence Hi n ^2 is 
empty). Moreover, given < t\ < ■ ■ ■ < t p < t and a pair 
(Qi, O2) in PJ -, we define the maps (C^ C T ) as above with 



t' — t 



^' — t 



and with t\ = t Pl , ■ ■ ■ , r- 



t 



We denote the set of all pairs of maps (C^,C T ) obtained 
in that way by T p j . 

With these definitions, we find the following exact for- 
mal expression for the reduced density matrix of the reg- 
ister 
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(l\p r (t)\m)=J2H) P dt 



p=0 



p 

o Jo 



tp />t2 

dtp— i ■ ■ ■ 



dti^i-iy E 

J=0 (Ci,C T )6T PJ 



(CoVoICi) ( II II(C T |A(4)IC T -i)(Ci-i|A(ti)ICi) ) cxp Li f [e( S ,C T M) - e( S ,C X M)] ds 

Vr=l S =l / L 



(11) 



^i[C] 



Here the elements of the reduced density matrix are ex- 
pressed in terms of the influence functional Z\ [(] , which 
describes the effect of the bath on the time evolution. In 
the next section we evaluate this influence functional. 



IV. THE INFLUENCE FUNCTIONAL 

To compute the influence functional Z\[Q we follow 
the method used for the spin-boson model in Ref. 22 . We 
define g\(s) — d x ll x (s) and denote by H x (s) the inte- 
grand in Eq. (9). Then, inverting d s and d x , we obtain 
ifiu(s) = d x {H x (s)U x (s)}, hence 

ig\(a) = H x (s)g x (s) + j( S )H re (C(s))Ux(s) • (12) 
Therefore, since g\ (0) = 0, we have that 

9x(s) = -iUx(s) f\l( S 'n( S ')H re (a s '))U x ( S ')d s ' . 
Jo 



As a consequence, defining 



(13) 



h x (s) = 



Tr, 



Z e Z x [(] 
we find that 

d X Z x {(] = 

hence 

Zi[C] =exp 



{U x (2t + (3)Ul(s)H rE (C( S ))U x ( S )} , 

(14) 
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j(s)h x (s)ds )Z X [(] 



1 r2t 

dX / j(s)h x (s)ds 
o Jo 



(15) 



(16) 



since by definition, Z [(] = 1. 

From Eq. (5) we can decompose h x as a sum by defin- 
ing 

/k+(s) = ^[cJ Tre { Ux{2t + ® u l( s ) b t u *( s )} > ( 1? ) 

and /k (s) by the same formula but with the annihilation 

operator 6k instead of the creation operator &£. Whence, 
we have 

hx(s) = ^2MC(s))fi(s)+MC(s)rf^(s). (18) 



Further, from the Heisenberg equation, defining & k (s) 
as [/|(s)6j c f/A(s), we obtain that 

iS s 6l( S )=7( S )^( S )[6L^ e + A^ e (C(s))]f/A(s). (19) 

Therefore, from the commutation relations for bosonic 
operators, we obtain for the following differential 
equation 



d.f£(s) - ±7(«) (^fe/k(«) - A0k T (CW)) 



(20) 



where (/> k+ (m) is given by Eq. (5) and </> k _ (m) is the 
complex conjugate of k+ (r7j). 

The solution of Eq. (20) is obviously given by 

/±(s)=e ±i ^W fe(0)+ 



T A^ S e^^( s ') 7 ( S ')0k T (C(s'))^' 



(21) 

Now, because of the invariance of the trace under per- 
mutations, we have the boundary condition f£(2t + fi) — 
/ k (0). Hence, from Eq. (21) we find that 



/k(0) 



Ae^/ 2 
2sinh(/3w fc /2) 



/' 

Jo 



[^(C^))-^^))] d«- 



Defining 



A k (s) = ^ k (C T (s))-^ k (C i W) 

s k (.s) = k (C T (s)) + k (C i W), 



(22) 



(23) 



where k is given by Eq. (5), we find that the influence 
functional is given by 



(24) 



with 

A b = Ecoth(/3w fc /2) 

f ds f ds'e- iWfc(s - s,) A k (s)*A k ( S ') 



k 

xRc 



(25) 



■5 



and 



I 6 =Vlm[/ da f ds'e- iw ^ s - s ')A k (s)*S k (s') 
k Uo Jo 

(26) 

Recall that and are constant on each interval 
[tj-i,tj[ with 1 < j < p + 1 (t = and t p+ i = t). We 
define 



4 = ^0, 



(27) 



and 

- to) 



to J to 



-1/ \ 1 
-(t-t ) + -5 



-iw fc (s-s') 



1 _ g-iwfc(t-to) 



(28) 



Note that 



/ ds / da' 

J tj —l ^ T7T. — 1 



(29) 



where 

-^fc™ = <7fc(*j _ *m-l) + Qk(tj-1 - t m ) 

~Qk{tj — t m ) — Qk{tj-1 — t m -i) 

Hence, we find that 

p+i 

A b = coth(/? Wfc /2) ^ Re - tj-i)^* A£ 



(30) 



p+i j-i 

EE Re 

j=2 m=l 



Af A™M^ m 



(31) 



and 



p+i 

* b =E£ Im hfe-^-i) A k*4~ 1 



j=2 m=l 



(32) 



V. CHARGE-QUBITS 



In order to evaluate the influence functional further we 
need to specify the form of the coupling in the Hamilto- 
nian (2). Therefore, we will focus on the example rep- 
resented in Fig. 1. We consider N equidistant qubits, 
with two electronic gates per qubit. We denote by q n 
the center of the qubit n (1 < n < N, we put qi =0). 
For qubit 1, q is the center of the upper dot and — q 



the center of the lower dot (q = |q |). Moreover, we 
write d = q 2 — qi, and d = |d|, the distance between 
qubits. Hence q n = (n — l)d, and q • d = 0. Each qubit 
has a single electron. 

For H re we first consider the coupling between elec- 
trons in the qubits and longitudinal phonons. This can 
be described by a Frohlich-type Hamiltonian 23 , where we 
have 



N 



-ik[(n-l)d+q aj] 



(33) 



71=1 



where cr™ is defined as <7™|Z) = l n \l)- This form describes 
the coupling g(k) between the charge of qubit n localized 
at (n— l)d+q cr™ and the phonon bath. Here, the spatial 
configuration is explicitly considered. 

Since the optical phonons are gapped at low frequency 
they do not contribute to low temperature decoherence, 
only acoustic phonons are relevant. We also consider 
only a linear coupling between the phonons and the 
qubits. Nonlinear couplings can also be included and 
it has been shown that these can be mapped to an ef- 
fective spin-boson model, albeit with decoherence rates 
(friction coefficients) that are very strongly temperature 
dependent 24 . Surface phonons are not considered either, 
since the quantum dots are usually embedded well in- 
side the semiconductor. Moreover, we assume that the 
phonons are only coupled to the diagonal operators of the 
qubit and we do not consider any non-diagonal couplings 
(involving terms like o x or <j y ) between the environment 
and the qubit. Such terms appear, e.g., when spin de- 
grees of freedom become relevant 2,21 . For charge qubits, 
nuclear spins are irrelevant and we restrict our study to 
the diagonal spin-boson model. In Sections VI, VII, and 
IX, we will further discuss the acoustic, electronic and 
deformation environmental degrees of freedom. 

Introducing the notations 



X(s) = 



[C T 00-C+(*)] 
[C T W + C X W] 



(34) 



we obtain that 

A k (s) = -2ig(k) sin(kqo) ]T e^"" 1 ^*) (35) 

and 

N 



N 



71=1 



£ k (s) - 2g(k) J2 e- ikd( "- 1} {cos(kq ) 



71=1 



-iX™(s)sin(kq )} 



(36) 

Now, as before, we denote £(tj_i) as f and x(t m ) as 
X™. Moreover, we denote by J the N x N Jordan bloc 
Jij = 6ji+i and write v = (l,e ikd , • • ■,e ikd ( lv - 1 '). With 
these notations, we find that 
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JV-l 



A£ A? = 4|. 9 (fc)| 2 sin(kq ) 2 (£'|f"> + E >e" ikdr + (£ m | J r |e> ikdr , 



(37) 



r=l 



and 



JV-l 



K S™- 1 = 4|. (fc)| 2 sin(kq ) 2 Ix™" 1 } + E ^ 1 ^lx m_1 >c- ikdr + (x m " VlfV 



kdr 



2|g(fc)| 2 sin(kqo)cos(kqo) ^ ikd/2 / _ ikdA A . • 
+ sin(kd/2) ^ 



Hence, introducing the following notations for r between and N — 1 and n between 1 and JV 

Q2±{t) = 2 E IsWI 2- 2 coth(/3^fe/2) sin(kq ) 2 [cos(kdr) — cos(u>kt ± kdr) =F o^/; sin(kdr)] 
k w fc 

Qi±(£) = 2 ^ |p(/c)| 2 — sin(kqo) 2 [sin^/j ± kdr) — w^t cos(kdr) =F sin(kdr)] , 



(38) 



*n(*) 

we find that 



2 V | 5 (fc)| 2 ^ [sin(u> fc i - kd(n - 1/2)) - sin(u, fe i - kd(n - 1/2 - N))] sin ( kq °) cos ( k( lo) 
k ^ 



(39) 



2 V \g(k)\ 2 — [cos(kd(n - 1/2)) - cos(kd(n - 1/2 - N))} 



sin(kd/2) 
sin(kqo) cos(kq ) 
sin(kd/2) 



JV-l 



Ecoth(/3o; fc /2)Re [a£* AJ?g fc (i)] = [WIOQS+W + (CI J r W)Ql-{t)] 

k r=0 

JV-l AT 

[Ai'sjr 1 ^)] - E ^ [(eiJ r ix m - 1 >gi + (t) + (x ro -V r iTOU*)] + E^ [*»(*) - *»(°) - m 



r=0 



n=l 



r 



_{40) 



where is defined as = 1 and v r = 2 for r > 1. In n 3 /L 3 ) with rij integers and |k| < a -1 , where a is the 
the next section, we show that in the continuum limit, lattice constant which of the order of a few angstroms. 
^f n (i) and $„ are zero and we provide compact formulas In the limit of an infinite volume, the sum over k can be 
for Q2±(t) and Q\±{t). Note that from Eq. (40), the replaced by and integral 
terms X b and A 6 of the influence functional (Eqs. (31) 
and (32)) can be expressed through the functions Qi±(t), 
Q r 2± (t), * n (t) and $„. 



E 



V S 



(27TC L ) : 



/•CL& 1 /"TT />27T 

/ duw 2 / d0sin(6>) / , 
Jo Jo Jo 



VI. LINEAR DISPERSION (ACOUSTIC 
PHONONS) 

In order to simplify the expressions for Q2±(t) and 
Qi±(t) we need to specify the dispersion relation. Since 
we arc interested in the effects of a phonon bath, the 
only relevant phonons at low temperatures are acoustic 
phonons, hence we assume a linear dispersion form ui^ ~ 
Ci|k|, where cl is the speed of sound in the sample. We 
consider a cubic sample of volume V$ = LiL 2 L 3 . The 
sum in H e and H re runs over all k = 2n(ni/Li, nil Li, 



(41) 

Further, we can assume that g(u>/ci,) decreases in such 
a way that in (41) we can replace the bound c^a -1 by 
+00. 

Recall that q and d are orthogonal. We choose 



kq = — q cos(0) 
cl 



kd 



CL 



-dsin(0) cos(0) . 



(42) 
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Since 



i' Tl 
Jo 



sin(y cos(0)) cos(y cos(6*))<y(sin(6>)) sin(0)cZ0 



y sin(yr) cos(yT)g(\/ 1 — r 2 )c?r = , 



(43) 

for any function g such that the integral exists, $„ and 
^ n (t) are equal to zero. 

We now compute the functions Q\ ± and Q\ ± . We first 
expand the expressions in cos(to k t ± kdr) and sin(u>fct ± 
kdr) in terms of cos(wfci) cos(kdr) and sin(wfci) sin(kdr). 
Further, because f(cos((f>))d(f) = for any odd func- 
tion /, the terms in sin(kdr) do not contribute to the 
integral. Therefore, we can omit the indices ±. We then 
introduce the transit time (r s ) and a dimensionless pa- 
rameter a 



d 2q 
t s = — and a = — — 

cl a 



(44) 



The parameter a represents the ratio of the size of a qubit 
over the distance between qubits (typically, a is smaller 
than unity). Integrating over <f) the term in cos(kdr) 
leads to a term 2ttJq (ujT s rsin(6)), where Jq denotes the 
Bessel function of the first kind. To perform the integral 
over 0, we use the formula (z > 0) 



sin (ycos(0))J o B (zsin(0))sin(0)d6> 



sin(z) sin (V z2 + V) 



\A 2 + 4y 2 



to find 
Q£(i) 



[1 - cos(wi)] coth{f3uj/2)duj 



Ql(t) = I [sinful) - wt] , 

Jo w 

where the spectral function is given by 
J r (ui) = auj 2 \g(uj I 'c L )\ 2 



(45) 



(46) 



sin(ur / aui q ) sin(uj\/r 2 + a 2 /aw q ) 
ur/aoj q u\Jr 2 + a 2 /au q 

~ (47) 

with ci = Vs/(2tt 2 Cl) and w q = (or,) 1 = c L /2q . The 
spectral function Jo(w) is defined by taking the limit r — > 
or 



J (u)) = ciu 2 \g(uj I 'c L )\' 



1 - 



sin(w/o> 9 ) 

Lu/Uj q 



(48) 



Our expression for J has the same form as the one ob- 
tained for a single double quantum dot 25 ' 26 . Moreover, 
apart from the term linear in t of Q\ (t) and the r depen- 
dence, Eq. (46) is identical to Eqs. (4.22a) and (4.22b) in 



Ref. 19 . More importantly, the particular r-dependence of 
the spectral function in Eq. (47) is the main reason why 
the decoherence induced by the bath of phonons does not 
lead to "superdecoherence" . Indeed, if we expand J r in 
terms of a/r we obtain 

J r (w) ~ -y (jy) \g(uj/c L )\ 2 cos(ruj/auj q ) , (49) 

when au>/(4ru> q ) <C 1. This ~ 1/r 2 dependence is re- 
sponsible for the suppression of "superdecoherence". In 
the unphysical limit, where a, the ratio of the size of the 
qubit and the distance between qubits, is much greater 
than one, there is no dependence on r for the expres- 
sions of Q\ and Q\ and the spectral function J r (u) equals 
Jo (a;), which leads to "superdecoherence". It is interest- 
ing to note that our result for the influence functional 
differs from that of references 5,6 because they have ne- 
glected the angle between k and d by introducing the 
"transit time" as ujt s — k • d. 

The influence functional can now be written in a com- 
pact form by first defining 



JV-l 



r=0 



Q b m (t)=J2»r(j r + J r )Q r m (t), 



(50) 



with m = 1, 2. Hence, 2Q^ n (t) are the diagonal elements 
of the N x N matrix Q b m (t), and 2Q r m (t), with r ^ 0, 
the off-diagonal terms. Following a similar notation as in 
Ref. 19 , we define 



a$* = QUtj - tfc-i) 



-Q 2 (tj 



- t k ) - Q\{t 



tk-i) j 



(51) 

and Xj k by the same formula, but with Q\ instead of Q\. 
Note that the part linear in t of Q\{t) does not contribute 

b 



to X h - U . With these definitions, we find that 



p+i p+ij-i 

a 6 = £<eiQ^ +EE^'i A ^ fe > 

i=i j=2 fc=i 

p+i p+ij— l 

x b = e™?& - ^-oix^h-EE^'i^i^- 1 ), 

j=l j=2 fc=l 

(52) 

and the influence functional is then simply given by Eq. 
(24), where A 6 represents the exponential decay due to 
the coupling to the bath and X b describes the phase. 

It is important to observe that for a single qubit (i.e. 
N = 1) we recover the usual formula for the influence 
functional of the spin-boson model with the spectral func- 
tion Jo(w). Expressions (52) are evaluated quantitatively 
in Sections VIII and IX. 



VII. THE COUPLING TO OTHER ELECTRONS 

The coupling due to the metallic gates can either be 
provided by the two dimensional electron gas (lateral 



8 



gates) or by metallic top gates. Each gate is consid- 
ered as a gas of free electrons. The gates are labelled by 
(n, u) with u = +1 for the gate above the qubit n, and 
u = — 1 for the gate below the qubit n (see Figure 1). 
The electron gas in the gate (n, u) is described as 



ktr J 



(53) 



ko- 



where /£" and /£" are fermionic operators. The gates 
are supposed to be isolated from each other, hence 
fermionic operators with different indices n or different 
indices u commute. We describe the coupling between 
the register and the electrons in the gates (n, u) as 



Tirf jj \ rnui rnu 

n nu ~ u nu / t Jko- Jk'<r 



(54) 



kk'cr 



where U nu is an operator acting on the register's Hilbcrt 
space 



N 



U nu =uJ2v(\3-n\)al 



(55) 



If we add the sum over n ranging from 1 to N and 
u = ±1 of H£ u + H^ u to the Hamiltonian (2), and com- 
pute again the reduced density matrix of the register, 
this leads to Eq. (11), where the bosonic influence func- 
tional is now multiplied by a fermionic influence func- 
tional Z'\Q which is a product of fermionic influence 
functional corresponding to each gate 



N 



= n n znu ^ 



(56) 



n— 1 u— ±1 



The electron-electron coupling is given by the Coulomb 
potential. Since each gate is much closer to the cor- 
responding qubit than the distance between qubits, we 
assume that V(r) = for r > 1 in Eq. (55). Hence, 
writing Vo for V(0), the interaction term reduces to 
U nu = uVo<Jz, an d we are left w ith a two- level system 
coupled to a fermionic bath by a contact potential. This 
model has been studied in Ref. 22 , where it was shown 
that for a density of states p(e) constant throughout the 
conduction band Ep of the bath, i.e. p(e) ~ p c~ € / EF 
(hence po = 1/Ep), the fermionic bath behaves as a 
bosonic environment with a spectral density of the ohmic 
form. This means, that aside from an adiabatic shift, 
i.e. a shift of the bias energy of the system, the influence 
functional Z nu [(] is identical to the influence functional 
of the bosonic bath, provided c/>k = .9fec" with gk real, and 
Sk dt^i 1 ^ ~ w fe) = ^(^0' w fth a spectral density function 
of the ohmic form 



J(oj) = rjuc 



-<»/»>i 



where 



• arctan 2 (7rp Vb) 



(57) 



(58) 



Here we use u>[ ~ Ep/h, where Ep is the Fermi energy 
of the bath. 

As a consequence, the influence functional becomes, 
Zf[Q = c- ix/ c- A/ , using Eqs. (23), (31) and (32), 
where X* and are given by Eq. (52) with the matri- 
ces Q{ 2 {t) = 8 id q{ 2 (t), where id denotes the identity 
N x N matrix, and q{ 2 (t) is defined as in Eq. (46) with 
the spectral function given in Eq. (57). The factor 8 
comes on one hand from the sum over u = ±1, and on 
the other hand from the factor \ in Eq. (34). 

VIII. DECOHERENCE FUNCTION 

We now compute the decreasing rate of the off-diagonal 
terms of the register's reduced density matrix p r {t). As in 
Refs. 5 and 6 , we consider the case where the dynamics of 
the register is trivial, i.e. A(t) = 0. This means that no 
quantum operation is performed. In this case, (l\p r (t)\m) 
is simply given by the term p = in Eq. (11) and can be 
computed exactly. Note that the extension to N qubits in 
operation (A ^ 0) is far from trivial 27 , and is discussed 
in Section X. If A(t) = 0, we find that 



(l\p r (t)\m) = (l\p r Q \m)e X p 



— i / [e(s,l) — e(s,m)]ds 
Jo 

xZ\[l,m]Z f [l,m] . 



(59) 



Moreover, we have 



1 



A 6 = ^{l - m\Q b 2 {t)\l - m) 

X b =- A {l-m)\Q\{t)\l + m) 
m\\ 2 q{{t) 



(60) 



A' =2\\l 
X* =0. 



For the diagonal terms (i.e. I = m), A 6 = and A-^ = 
0, whereas for the off-diagonal terms, with Eq. (50) we 
have 



A" 



1 



Ql{t)\\l-mf 



^ (l-m\J r \l 

r=l 



m) Q5(t) 



- TO 



Ql{t) 



(61) 



Therefore, for the most off-diagonal terms (i.e. I 
2(1, • • • , 1)), we find 

A 6 =2NQl(t)(l + e(t,N)) 

X b =0 

Af =8Nq f 2 (t), 



m 



(62) 



with 



JV-l 



e(*,A0 = 2£(l--l) 



Ql(t) 



(63) 
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which is small and bounded in the variable N since Q r 2 ~ 
r~ 2 as discussed in section VI. This implies that the 
decreasing rate is proportional to the number N of qubits 
in the register. Note that if we neglect the r dependence 
and put Q\{t) = q 2 (t), for all r > 1, we find the result 
stated in Refs. 5 ' 6 . 

The decoherence functions evaluated above are the 
ones corresponding to the most off-diagonal elements of 
the density matrix. They correspond to the maximum 
decoherence rate of the register. For other matrix ele- 
ments, the bound |(£|J n |OI < U\\ 2 ieads to 



&-(*)< \(l\ P r (t)\m)\<b + (t), 



(64) 



with 
b±(t) 

and 



m) | exp 



-Q 2 (t)\\l- m \\ 2 (l T e(t,N)) 



x exp 



-211/ 



|2„/ 



(65) 
(66) 



Note that e(t, N) > for all t and N. If there is e < 1 
such that e(t, N) < e for all t and N, the expressions and 
bounds in Eqs. (64)-(66) imply that there is no superde- 
coherence nor a decoherence-free subspace. This result 
is consistent with the additivity of decoherence measures 
in the short time limit. 28 



IX. PIEZO AND DEFORMATION PHONONS 

The coupling function g(k) in Eq. (33) depends on 
the nature of the phonon coupling. In this section, we 
discuss two important cases: the phonon interaction cor- 
responding to the deformation potential and the piezo- 
electric phonon interaction. 

In experiments on single GaAs/AlGaAs double quan- 
tum dots it was argued that the main contribution to de- 
phasing is due to the piezoelectric phonon interaction. 29 
In this case it was shown that for a double quantum dot 
the spectral function is given by 25 



J(u) = guj 



1 



sin(w/o; g ) 



Uj/uO q 



ij/u 



(67) 



where e~ u / u * is the high frequency cut-off function, 
which could have different forms. Hence, comparing Eqs. 
(67) and (48) we find that the piezo case corresponds to 
taking c\\g(uo / cl)\ 2 — $-e~ w / u c. In recent experiments 
on double quantum dots it was found that g ~ 0.03 7 ' 29,30 . 

For the deformation potential the dependence is given 
by 25 



J(u) 



UJ 



1 



sin(u;/u; g ) 

Uj/Ujq 



(68) 



where w 2 ~ 10 25 s 2 for typical GaAs values 26,31 . This 
is exactly the form we obtain by taking c\\g{u) / 'cl)\ 2 — 



in Eq. (48). In the limit 



0, both 



behaviors can be related to the widely used parameter- 
ized form of the spectral function J(u>) ~ uj s e^ u / UJ <^. 
Therefore, the ohmic case s — 1 is the analogue to 
the piezoelectric phonon interaction and the superohmic 
case s = 3 is the analogue to the deformation potential 
phonon interaction. 




FIG. 2: Maximum decay functions for N = 10 qubits given 
by Eqs. (62) at zero temperature for piezo and deformation 
phonons and electronic baths. The constants are g = 0.03, 



= lO^s-^, 
400nm, uj b r 



g = 9.3 • 10" 8 , c L = 5 • lO^m/s, g =50nm, 
= uj a — 5 • 10 10 s _1 (phonons), and lo[ = 1.3 ■ 



10 15 s 1 (electronic bath). The corresponding dotted lines are 
2NQ2' ptez ° and 2NQ^' de ' f ' , which illustrate the expressions 
for a single qubit. 



We evaluate A fc , A' and Q r 2 using Eqs. (62), (46)- (47), 
(67)-(68), and (57) for typical values of a multiple coupled 
quantum dots system imbedded in GaAs/AlGaAs for the 
case where we have N = 10 3 qubits. The results are 
shown in Fig. 2. The cut-off frequency is important, since 
it defines a characteristic time scale. For the electron- 
phonon coupling the relevant phonon frequency is given 
by the smallest extent of the electronic wave function in 
the quantum dot, which we assume to be Ao=100 nm. 
Hence, oj h c = cl/Xq- This also implies that the coherence 
time decreases with a stronger quantum dot confinement. 
In contrast, lo[ for the electronic bath coupling is given 
by the Fermi energy of the gates. For lateral gates the 
Fermi energy is given by the two-dimensional electron 
system, where we assume a typical Fermi energy of 8.9 
meV, which corresponds to a Fermi wavelength of 50nm. 
The coupling constant is given by (58), where Vq is the 
charging energy of the dot, which we assume to be 1.2 
meV and corresponds to a typical dot to gate separation 
of lOOnm. This leads to rj = 3.2 • 10~ 2 . For a similar 
geometry but with top metallic gates, the Fermi energy is 
about 5.5 eV (for gold) and the charging energy is similar 
to the lateral gates geometry, i.e., 1.2 meV. This leads to 
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r] = 9.3 • 10 -8 . The values of go and d are relevant to 
recent experiments 29 . With these parameters we obtain 
a decoherence function for the lateral gates (not shown 
in Fig. 2), which is five orders of magnitude larger than 
when we consider only top metallic gates (shown as A-^ 
in Fig. 2). Hence, we will only consider the top gate 
geometry in the remainder of this discussion. 

The main contribution to decoherence is clearly given 
by the piezo phonons as was argued earlier 29 . The 
coupling to the electronic leads (metallic) introduces a 
smaller decoherence decay and the form of its time- 
dependence is the same as the single qubit case (except 
for the prefactor). For the phonon bath, considering N 
qubits, instead of one qubit, modifies the form of the time 
dependence in addition to the prefactor. The difference 
in behavior is illustrated by the solid and dotted lines in 
Fig. 2. 

The small oscillations seen in the figure are reminis- 
cent of coherence revival 32 and are most likely due to 
phase exchange between the qubits via the environmen- 
tal bath. The saturation of A 6 at large times is similar to 
the saturation seen in the superohmic case of the spin bo- 
son model. 19 This saturation occurs because of the small 
density of low frequency modes in the spectral function. 
Indeed, at low frequencies, the leading order of the spec- 
tral function for piezo phonons in Eq. (67) is given by 
J(u>) ~ uj 3 (superohmic at low frequencies). Eventually, 
full decoherence would occur if we include the exchange 
of energy between the qubits and the bath (A ^ 0). This 
introduces another time scale T± above which all coher- 
ence is lost. However T\ is usually much longer than 

u) b c , the typical time-scale of the decoherence due to 
quantum fluctuations. We leave the discussion of the en- 
ergy transfer processes to Section X, where we consider 
A ^ 0. 

The decoherence time due to the quantum and thermal 
fluctuations (non-dissipative) can be obtained from Eq. 
(62) and is given by A b (tdec) — 1. Hence, from Fig. 2 we 
can estimate that tdec — 5ps for N = 10 3 qubits at zero 
temperature. The temperature dependence of the decay 
function A & (Eq. (62)) for the coupling to the dominant 
piezo phonons is shown in Fig. 3 for t=l, 10, and lOOps. 
As expected, the main effect of an increasing tempera- 
ture is to increase the decay function. At low tempera- 
tures the decay function saturates close to lOOmK and 
the decoherence mechanism is only due to quantum fluc- 
tuations. It is interesting to note that for a small number 
of qubits, for example N = 10, the decoherence function 
at Ins is only 10%, which means that there is very little 
decoherence. However, this small decoherence would still 
lead to an error rate during a quantum operation. 26 

The temperature dependence of the decoherence func- 
tion for a single qubit, which is described by the function 
Q0, P iezo^ j g g j 10wn m foe mse t f pig 3 m dotted lines. 

This shows that the overall dependence is similar at low 
temperatures and small times, where 2NQ 2 is actually a 
good approximation to A b . For higher temperatures and 
longer times the exact expressions (62)-(63) have to be 




t[s] 

FIG. 3: Time dependence of A b piezo (Eq. (62)) for piezo 
phonons and N = 10 3 qubits for different values of the tem- 
perature. Inset: Temperature dependence of Ap iezo for 1, 10 
and lOOps. The corresponding dotted lines are the functions 
2iVQ°' pl "°. 

used. In Fig. 4 we plot the functions Q\ and Q 2 from 
Eq. (46) for the first values of r. Hence, in practice, it is 
enough to sum over a few values of r in order to calculate 
A b . For the phase function X b and Q\ the situation is 
very similar. 

t[s] 

10" ,a 10"" 10"'° 10"' 
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4NQ 2 2 /~^; 
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/4NQ 2 3 















10* 
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-4NQ t ' 


-4NQ, 2 
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io" io"" io"'° 10" 
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FIG. 4: The functions 2NQ° 1 , 2NQ° 2l ANQi and ANQ2 for 
piezo phonons and r — 1, 2, 3 and N = 10 3 . 



X. DYNAMICS IN THE REGISTER 

In previous sections VIII and IX the results were ob- 
tained assuming that no quantum operations are per- 
formed (i.e. A = 0). In this section we introduce non- 
trivial dynamics in order to study the effects of quantum 
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operations on the decoherence rates. 

We assume that A is constant. Suppose that for a 
set S of qubits, the dynamic is trivial (i.e. the quantum 
operation does not involve those qubits). We take (£,x) 
to be a given path occurring in Eq. (11). Define £* r as 
£* r = if n is not in 5, and £^ = £ n (0) otherwise, and 
£ dy as £ — £* r . Define x tr and x dy m the same way. The 
subscript tr stands for "trivial" , and dy for "dynamical" . 
Then, we find that the influence functionals for bosons 
and fermions depending only on the trivial part of the 
path are given by 

A t y ={e\Q b 2 f (t)\e r ) 

(69) 

x b j =(e\Q\ J (t)\x tr ). 

Note that since Q{(t) is a diagonal matrix, X( r = 0. 

We now define q\'2 as the off-diagonal part of Q^. 
By definition, we have (£, tr \£ dy ) = 0. As a consequence, 
we find that the part of X b depending on cross terms 
between the trivial and the dynamical parts of the path 
is given by 

x b Jj2{e vj \Q\{t,)-Q\{t^)\x tr ) {7Q) 

+(e r \Ql(t-t j _ 1 )-4'i(t-t j )\x dyj - 1 ), 

whereas x* = since Q{ is a diagonal matrix. The cor- 
responding term for A 6 , which we denote by X b , is given 
by the same formula with x* r , x 1 an d i replaced 
by £* r > £ j and 2 respectively. For the same reason as 
above, A-' = 0. We are interested in the dynamics of the 
register's reduced density matrix for time scales smaller 
than the decoherence time, which we have estimated in 
the previous section (t < tdec)- As a consequence, since 
for t < tdec, Qi(t) and Q%{t) are essentially zero for all 
r > 1 (see Fig. 4), we can neglect the cross terms of the 
bosonic influence functional, that is we can safely assume 
that x b = and A b = 0. 

In general, a quantum computation can be achieved 
by successive single-qubit and C-NOT operations. For 



C-NOT gates, we consider the set-up proposed in Ref. 8 , 
described by the Hamiltonian 

-A — - — <g>a x -e — - — ®(T Z , (71) 

where a NOT operation is achieved on the right qubit af- 
ter a period given by AT no t — § , whenever the left qubit 
is in the state | — 1). Note that in a semiconductor charge 
quantum register, only C-NOT operations between near- 
est qubits can be achieved. For single-qubit gates, we 
consider the Hamiltonian — Aa x ~ ea z . Three single-qubit 
operations are needed for a quantum computation. Two 
of them are trivial, i.e. A = and T\s — 3tt/4 and 
T 2 e = 7n/8. The last one is the Hadamar gate given by 
A = e and T Had VA 2 + e 2 = § . 

As a consequence, the decay rate of the most off- 
diagonal terms of the register's density matrix for a single 
C-NOT operation or qubit rotation, will remain propor- 
tional to N. For parallel qubit rotations or C-NOT op- 
erations, the compact formulas for the bosonic influence 
functional are a good starting point to discuss a gener- 
alized NIBA approximation. However, a more detailed 
analysis is beyond the scope of this article. 



XI. CONCLUSION 

We have analyzed the decoherence process in a solid 
state quantum register with N qubits. We showed that 
the decay rate of the most off-diagonal terms of the regis- 
ter's density matrix is proportional to N in all situations 
relevant to a scaled charge solid-state quantum computer, 
where the qubits are coupled to a common phonon bath 
and to independent electronic gates. We obtained com- 
pact expressions for the iV-qubit decoherence functions 
and argued that when performing quantum operations 
the decoherence function follows a very similar depen- 
dence as compared to the static case. 
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